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We investigate the steady state photon transport in a nonequilibrium collective-qubit model. By 
adopting the noninteracting blip approximation, which is applicable in the strong photon-qubit 
coupling regime, we describe the essential contribution of indirect qubit-qubit interaction to the 
population distribution, mediated by the photonic baths. The linear relations of both the optimal 
flux and noise power with the qubits system size are obtained. Moreover, the inversed power-law 
style for the finite-size scaling of the optimal photon-qubit coupling strength is exhibited, which is 
proposed to be universal. 

PACS numbers: 05.60.Gg, 42.50.Nn, 42.50.Lc 


I. INTRODUCTION 

Deep understanding and optimal control of quantum propagation in low-dimensional light-qubit (atom) hybrid 
systems is of fundamental interest and practical importance [1]. Particularly, the collective effect on the information 
and energy transport due to light-qubit (atom) scattering, as a novel measuring feature, has recently attracted dramatic 
attention M. Many works have been carried out to observe such effect, ranging from solid state physics [6, 7], 
quantum biology [8]. to quantum optics [9, So]. The transport scheme, quantum flow from hot source to cold drain, 
can be typically established by applying the thermodynamic (e.g., temperature) bias, in accordance with the second 
law of thermodynamics. 

The prototype paradigm of collective-qubit systems to characterize quantum transport based on photon-qubit 
interaction, is termed as Dicke model. It was originally pointed out by R. H. Dicke, which described N identical 
two-level atoms coupled to single radiation field mode [llj. Dicke model has been extensively studied in quantum 
superradiant phase transition associated with collective self-organization of atoms, which predicts the universal finite- 
size scaling effect 13 El. Moreover, the influence of photon dissipation on the collective-qubit model is analyzed by 
considering continuous radiation modes. The anomalous superradiant-like relaxation, dynamical quantum beat and 
quantum phase transition have been unraveled, which significantly differ from the counterpart in spin-boson model 

(AT = 1) PIl. 

Recently, steady state transport behavior of collective qubits weakly coupled to two photonic baths is studied 
together within Redfield scheme and full counting statistics [13, H8| , where the thermodynamic bias is applied to drive 
the unidirectional photonic flow ESE1- The influence of bath temperatures on the generation of collective quantum 
transport is intensively analyzed. As is known in nonequilibrium spin-boson model, strong system-bath coupling 
plays nontrivial role to exhibit nonmonotonic flux behaviors, by including the nonadditive and nonresonant transport 
processes between qubits and bosonic (e.g., phonon) baths |21 25j|. Hence, it is indeed desirable to study the effect of 
strong photon-qubit interaction on quantum transport in collective-qubit system. 

In the present work, we adopt the nonequilibrium noninteracting-blip approximation (NIBA) |2j], HI, HHHl to 
study the steady state collective transport behaviors of multi-qubits system in strong coupling regime. The indirect 
qubit-qubit interaction mediated by photonic baths, is clearly revealed, contributing to collective behaviors. Then, the 
effect of strong photon-qubit coupling on photonic energy current fluctuations (e.g., flux, noise power) is investigated, 
and the finite-size scaling of the corresponding optimal observables is discussed. This paper is organized as follows: 
in section II, we describe the collective-qubit model and derive the quantum kinetic equation in collective-angular 
momentum basis. The comparison of steady state population in strong coupling regime is made with the counterpart 
within Redfield scheme. In section III, the full counting statistics of photonic energy flux is established, and the 
finite-size scaling of the current fluctuations is analyzed. A concise summary is given in the final section. 
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FIG. 1: (Color online) Schematic description of nonequilibrium qubits system: the left red arch and right blue arch represent 
the source and drain photonic baths, with the temperature given by Ts and Td, respectively; the central pink circles combined 
with the arrows describe the two-level qubits; the cured red and blue lines describes the system-bath interaction. 


II. NONEQUILIBRIUM COLLECTIVE-QUBIT SYSTEM 

A. Model 


The nonequilibrium collective-qubit model at Fig.[lJ composed by N identical two-level qubits interacting with two 
photonic baths, is described by 

N . N N 

^ = -§E^ + yE^ + EE K, k ^(bi, v +kv )+E (i) 

n= 1 n= 1 n= 1 v,k k,v 

where a™ and a™ are the Pauli matrix representing the nth two-level qubit, eo is the Zeeman splitting energy of the 
qubit, A denotes the coherent tunneling strength between two levels of the qubit. b\ v (bk,v) creates (annihilates) 
one photon with the frequency LJk and momentum k in the vth photonic bath, and A n ,k,v describes the coupling 
strength between nth qubit and photon with momentum k in rdh bath. In this paper, we simplify the system-bath 
coupling independent of the specific qubit ^ kv = ^k,v The interaction of the qubits system with the photonic baths 
is characterized by the spectral density J v (uj) = 47r^ fc |Afc jV | 2 5(uj — Uk)- Here, we assume that the typical Ohmic 
function characterizes baths p9l-[3ll|. specified as J v (uj ) = a v uoe~ u / Uc ' v , where a v is the coupling between the qubits 
and vth bath, and cj c?v is the cutoff frequency of the vth bath. 

For the qubits do not interact with each other, we introduce the angular momentum operators J z = \ 1 °z 

and J x = \ J2n=i to describe the collective behaviors of qubit system. Then, the nonequilibrium collective-qubit 
model is simplified to 

H = e 0 Jz T AJa, + 2 J z A v,k(Pk,v ^k,v) T ^ ^ ^k^\ v ^k,vy (2) 

v,k k,v 


where the angular momentum operator J x = d(J + -f J_), with J± creating and annihilating angular momentum. 
They are described as J±\j,m) = V~gm\j,rn) under the angular momentum basis \j,m), with j = N/2 and = 
j(j + 1) — m(m±l). They obey such commutating relations [J + , J_] = 2 J z , and [J±, J z \ = ±J±. 

As two photonic baths have the temperature bias (e.g., Ts > Td), there exists unidirectional photonic flux from hot 
source to cold drain. Here, we are interested in the influence of the multi-photon excitations on quantum transport, 
which usually occurs beyond weak system-bath coupling regime, it is practically useful to transform the model at 
Eq. ([2]) by using the canonical transformation to Ht = P^ HP, with the unitary operator P = e lJzB and the collective 
photonic momentum operator 


B — ^2 B v - i V 


2 (i t 




(bk.v bk, v )- 


(3) 


v=S,D k,v 

The transformed Hamiltonian is shown as H t = H s + V s b + v 0Jkb\ v bk,v The transformed system Hamiltonian is 


H-s — e 0 Jz 


( 4 ) 
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A 2 

where £ = 4 v is the reorganized energy, contributing to the collective excitations of qubits system. Given the 
Ohmic spectrum, the reorganized energy is specified as 


£ = 



( 5 ) 


which is linearly proportional to the system-bath coupling strength. The eigenbasis is | j, m) for the Hamiltonian at 
Eq. (J4]) , and the eigenvalue is expressed as H s \j,m) = Ey m |j,ra), with £/ ?m = —eo m — £ra 2 , m = —j, —j + 1,..., j. It 
is interesting to note that for the single qubit case (.N = 1), the collective-qubit model is reduced to the well-known 
nonequilibrium spin-boson model [UMil- The photonic bath mediated long-range coupling term — £ J 2 is simplified to 
the constant parameter—£/4, which can be safely ignored to study the transient dynamics and steady state behaviors. 
For two or more qubits ensembles, the influence of this collective term on the dynamics becomes apparent beyond the 
weak system-bath coupling regime, mainly due to the reorganization of the system energy levels Ej :Tn . Whereas such 
feature can not be unraveled from the spin-boson model. 

The transformed system-bath interaction is given by 

Vsb = 2 ^J+ + e 2 ^J_). (6) 

It is clearly shown that the photon transfer process is contributed by the excitation of the collective-qubit system 
from lower state to higher one with absorbing photons from the baths, or relaxation of the qubits from higher state 
to the lower one by emitting photons to the baths. Moreover, spin flips are accompanied by multi-photons collective 
transport. This intrinsic picture may be difficult to uncover from the original Hamiltonian at Eq. 0 . 


B. Quantum kinetic equation 


We apply the nonequilibrium noninteracting blip approximation (NIBA) (HI [HI HH-Hil to study the collective- 
qubit system dynamics. Born approximation is included to perturb the transformed Hamiltonian at Eq. ([6]) up to the 
second order. It is known that this approximation is valid for weak tunneling A<cj c ^, at high temperature and/or 
strong system-bath coupling regime [2l|, where cv C;V is the cutoff frequency of the vth bath. Moreover, we consider 
the Markov approximation, based on the assumption that the relaxation time of photonic baths is much shorter than 
the counterpart of the qubits system. Then, the quantum kinetic equation is expressed as 


^ — (Am—1 "b ft/j-Pm + ^m-l^m-1 + 1? (7) 

where P m = (j,m\p s (t)\j,m) is the qubits population at angular momentum state | j,m). The transition rates 
excite (relax) qubits from state \j,m) to |j ,m + 1) (from state |j, m + 1) to |j,m)), demonstrating the cooperative 
transfer process between the qubits and photonic baths. They are expressed as 


K 


± 

m 



duC±(w)du, 


( 8 ) 


where C±(cv) = Cs(Toj)CD(Eoj^A m ) is the transition kernel, and the energy gap between level m and m + 1 is 
A m = P m +i — E m . C v {uj) = J^° oo dte luJt ~ ( ^ v ^ (v = S,D) is the probability density of the vth bath, describing 
the absorption of energy uj from the vthe photon bath if a; > 0, or the release of energy — u if uj < 0. Q v (t) = 
) 2 [{^ n k,v + 1)(1 — coster) + isinej/cT] is the photon propagator from vth bath, which is obtained by the 

thermodynamic statistics of collective photon momentums — ( e -*B v (t) e tB v (o)y found that the local vth 

transition density fulfill the detailed balance relation as C v (—uj)/C v (uj) = e~ f3vUJ . However, the transition kernel 
C ± (cv) breaks this relation, due to the nonequilibrium bias (Ts^Td)- This implies the complex transport processes, 
cooperatively contributed by two photon baths. Moreover, it is found from Eq. m that the transition process is 
described by the perturbation to the lowest order of the coherent tunneling A 2 , whereas photons are non-perturbatively 
involved in the energy transport process. 

Since we mainly focus on the physical behaviors in the high temperature T$(d) > e o an d/or strong system-bath 
coupling regime. The well-known Marcus limit is considered for calculating simplification (32|. Briefly, it can be 
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FIG. 2: (Color online) Comparisons of steady state population distribution ((a) and (b)) and eigen-energy levels ((c) and 
(d)) within nonequilibrium NIB A at Eq. 0 and Redfield scheme at Ref. 0. The parameters are given by N = 6, eo = 0, 
(jj c — 10A, a — 0.1A, Ts — 4A and Td — 2A. 


achieved by expanding the photon propagator Q v (r) in short-time regime as Q v (r) 
coupling strength £ v = a v uj c , v /n. Consequently, we obtain the transition kernel C v {u) 
the transition rate 


K 


± 

m 


(-) 2 a+ / 1 r (A ro ±&±fr?) 2 1 

1 2 > 9m \j Tsts + T d ^d P [ 4 (Tstis + T D £ D y- 


£, v (t 2 //3 v + ir), with the 
f ^ c (& £v) and 


(9) 


By inserting this expression of rates into the equation of motion at Eq. we are able to analyze the collective steady 
state behaviors with low computational cost. 

Unlike the Redfield scheme in the weak coupling regime, which shows resonant and additive transport behaviors IS 
[ 20 ], population dynamics at Eq. © describes the nonresonant photonic energy transport process [HI, H3, HH, which 
is clearly demonstrated from the transition rate at Eq. ©• Specifically, demonstrates that as the collective-qubit 
system absorbs energy A m through the transition from the state | j, mn) to |j, ra+1), the source bath will contribute uj, 
and the drain bath will supply the rest A m — u. While means that as the qubits release energy A m from the state 
\j, m + 1) to | j, m), the source bath will gain the uj and the drain bath will absorb the rest A m — uj. Moreover, it should 
be noted that a similar derivation of the quantum kinetic equation for the multi-level model was obtained [HI , which 
focused on the quantum fluctuation theorem. In this paper, we mainly exploit the influence of the strong system-bath 
coupling on collective transport of the finite-size nonequilibrium qubits system. 


C. Steady state population 

We study population distributions after long time evolution lirn^oo dP< 2^ = 0. From the equation of motion for 
population dynamics at Eq. 0, the population relation of nearest-neighboring states is given by P m+ i/P m = 

For the single photonic bath, it is reduced to P m+ i/P m = e -/3Am , due to the detailed balanced relation of the transition 
processes. While under the thermodynamic bias (i.e. Ts > To), such local balance will break. Unlike the Redfield 
scheme, where the effective ’equilibrium’ state can be quantified through the additive feature of photonic baths 0, 
it is difficult to analytically obtain such effective distribution in strong coupling within NIBA, due to the nonadditive 
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transfer behavior of photonic baths. However, for small reorganized energy £<Ceo at off-resonant case, the energy gap 
can be approximated by A m = eo + (2 m + l)£~eo, which becomes level independent. Hence, the transition rates 
can be simplified to y±, with the rate kernel y± = (A/2) 2 ^ f^° oQ dujC±(uj) and the probability density 

C±(uS) = (7 i s'(=FCo’)C j d(±co’=fA). To show the expression of population distribution, we assume P m = ay m , and replace 
it to the kinetic equation at Eq. resulting in y = 7 + /y_. Moreover, through the conservation of the probability 

P™, = 1- the constant parameter is obtained as a = i with j = N/ 2. Hence, the population distribution is 

shown as 


P — 

J- m. — 


y N / 2 {\-y) 

l-y N +l » ■ 


( 10 ) 


It should be noted that though the profile of the steady state population at Eq. (ITU]) within nonequilibrium NIBA 
seems the same as in the Redfield scheme jl9j . the transition ratio y from two schemes are quite different, which 
results from the distinct transfer pictures [24| . 

For strong system-bath coupling, by applying the nonequilibrium NIBA approach we numerically plot the steady 
state population at resonance case (eo = 0), as shown at Fig. |2ja). It is found that population shows non-monotonic 
behavior in eigenbasis |j, m), with maximal probabilities at two highest excited states \ j, ± j). Moreover, it is symmet¬ 
rically arranged from m = 0. In sharp contrast, the population distribution within the Redfield method is expressed 
as the monotonic decrease at Fig.[2fb), with the maximal probability expected at ground state |j, — j). The difference 
of population based on two schemes mainly originates from the long-range interaction of qubits, mediated by photonic 
baths. Specifically, at strong coupling regime, the influence of many photon excitations on the indirect qubit-qubit 
interaction becomes apparent. This results in the nonlinear term — at Eq. which is unraveled by the nonequi¬ 
librium NIBA. Then, the energy levels of qubits system are changed to E m = —£ra 2 shown at Fig. |]Jc), with lowest 
levels at states |j, ±j). Considering the steady state transition balanced relation ft+P m = /€“P m +i derived from 
Eq. ©, and the expressions of transition rates at Eq. (|9]), we are able to obtain that two highest excited state occupy 
the maximal probability and state \j, 0) has the minimal probability. However, the energy levels in Redfield scheme 
Em = Am increase with angular state m at Fig. |2fd), which means that the steady state population should decrease 
accordingly. Hence, we conclude that the indirectly collective interaction between qubits is crucial for steady state 
behavior in strong system-coupling regime. Whereas it can not be exploited by the Redfield scheme, which is usually 
applicable in the weak coupling limit. 


III. FULL COUNTING STATISTICS OF PHOTONIC ENERGY FLUX 


By applying the full counting statistics HzUH, we next obtain the expression of cumulant generating function 
of the nonequilibrium collective-qubit system, to count the photon energy transport from hot source to cold drain. 
By including the counting field parameter to count photon energy into the drain bath, the Hamiltonian at Eq. m 
is changed to H(x) = e lxHD ^ 2 He~ lxHD ^ 2 ^ with Hu = the drain bath Hamiltonian. Then, under 

the modified canonical transformation, with the counting parameter embedded unitary operator P x = e lJzBx and 
the modified collective bath momentum B x = i^ kv 2X J k v e lxuJkSv ’ D / 2 b k v + P.c., the Hamiltonian is modified to 

H t (x) = H s + Hb + V 3 b(x)- The reorganized system-bath interaction is described as V s b(x) = y (e~ zBx J + + e zBx J-). 
Similar to the procedures to obtain the Eq. o under the Markovian-NIBA approximation, the quantum kinetic 
equation combined with the counting field parameter is given by 


dP& 

dt 


(Pm—1 "h ^m)-^ > m "b ^ rn —\ (yX)Pm—l "b 


where the modified transition rates are 


(ii) 


A f‘°° 

4(x) = (-) 2 gle^ A "**y 


( 12 ) 


In absence of the counting field parameter, P x ~° are the state populations P m , and Eq. (fill) is reduced to the standard 
quantum kinetic equation at Eq. ©. Meanwhile, the modified transition rates return back to Eq. (j8j). In the Marcus 
limit, these transition rates at Eq. (IT2l) are simplified to with shown at Eq. m and the 

counting factor 


F mU) = exp{=F*A m x - 


TsTpljsiD r . , 1 , 
Tsis + Toin [ X( T S + 


TA rn - tD 
Td£d 


) + x 2 ]}- 


(13) 
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FIG. 3: (Color online) (a) Photonic energy flux of finite-size nonequilibrium qubits system in wide system-bath coupling regime; 
(b) finite-size scaling of the optimal energy flux: J opt — max a {J}; (c) finite-size scaling of the optimal system-bath coupling, 
which corresponds to the optimal flux. The parameters are given by eo = 0, uj c — 10A, Ts — 4A and Td — 2A. 


Then by re-arranging populations in column style |P(x, t)) = [P* ■, , • ••, P*] T , Eq. (JTTJ) can be re-expressed as 

J^-|P(x, t)) = —£(x)|P(x,t)), with £(x) tri-diagonal matrix dependent on the transition rates. Hence, the generating 
function can be formally obtained as G(x) = lim*-^ \ In (/|P(x> t)), with the ^ un ^ vector (I | = [ 1 , 1 , 1 ]. For 
the long-time limit behaviors are of primer interest in this study, the corresponding cumulant generating function is 
practically approaching the ground state energy Eq(x) °f transfer matrix £(x), having the smallest real part. Based 
on the cumulant generating function, the nth order cumulant of photonic energy current fluctuations can be steadily 
obtained as lim^oo (( Q n ))/t = lx=o? and the steady state energy flux is the lowest order J = lx=o- 


A. Steady state energy flux 

We study the influence of system-bath coupling on the steady state photonic energy flux beyond the weak interaction 
regime, as shown at Fig.(3ja). By increasing coupling strength, the energy flux is firstly raised to reach the summit. It 
mainly results from the contribution of the nonresonant photon transfer processes. If we further enlarge the interaction 
strength to the ultra-strong coupling regime, energy flux is suppressed due to the detrimental photon-qubit scattering. 
As a result, the non-monotonic feature appears, which provides guide to the optimal control of the photonic flux. 
This optimal behavior is similar to the nonequilibrium spin-boson model within NIBA scheme Hill , which is the 
special case of the collective-qubit model by setting N = 1 . Next, we analyze the finite-size effect of the energy flux. 
By increasing the system size, the optimal value of the energy flux is accordingly amplified. Specifically, the optimal 
flux shows linear scaling with the system size J opt (xN , shown at Fig. (3^b). Moreover, the optimal coupling strength 
shrinks with the increase of the system size at Fig. m c), exhibiting inversed power-law scaling a opt ocN 7 , with the 
scaling exponent 7 = 2.00T0.06. It should be noted that these novel features can not be observed within the Redfield 
scheme jl9|, where the resonant energy flux linearly increases by enlarging the system-bath coupling, and no optimal 
behavior emerges. 

Then, we turn to study the finite-size scaling effect of the energy flux under high temperature bias (e.g., Ts = 
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FIG. 4: (Color online) (a) Photonic energy flux of finite-size nonequilibrium qubits system with different thermodynamic biases; 
(b) finite size scaling of the optimal flux: J opt — max a { J} with bath temperatures Ts — 20 A and Td — 2 A; (c) finite size scaling 
of the optimal system-bath coupling with bath temperatures Ts — 20 A and Td — 2 A, which corresponds to the optimal flux. 
The other parameters are given by eo = 0 and uj c — 10A. 


20 A, Td = 2 A) to find the universal feature, shown at Fig. |U It was reported that within Redfield weak-coupling 
limit, the energy flux in the high temperature regime asymptotically reaches summit, and superradiant transport 
signal is found(Joe N : 2 ) M- However, in strong coupling regime, the energy flux exhibits nonmonotonic behavior, 
which demonstrates disappearance of the peak platform in high temperature bias. Besides, its optimal value shows 
the linear scaling with the system size, which implies that the linear form may be universal for the optimal photonic 
flux. This results seem to be a sharp contrast to the counterpart in the weak coupling regime. However, the deviations 
in the same quantum system between two schemes just reflect the different physical manifestations. Hence, they are 
complimentary to each other to fully unravel novel system features. Moreover, the coupling strength corresponding 
to the optimal energy flux shows inversed power-law case <a o:P £ 0 c 7 V _7 , 7 = 2.00T0.10. This strengthens the proposal 
that the finite-size scaling of the optimal coupling strength may be universal. 


B. Noise power 

We study the influence of system-bath interaction on the noise power, defined by S = d 2 G(x)/d(ix) 2 \ x =o- By 
tuning the coupling strength into strong regime, the monotonic behavior of the noise power is clearly exhibited given 
by Fig. 01 a), which is similar to the photonic energy flux at Fig. [3ja). To further describe the statistic feature of 
fluctuation of the energy flux, we introduce FF = (( Q 2 ) — ( Q) 2 )/(Q ) = S/J with Q = cdkb\ d H,d^ defined as the 

ratio of the standard deviation to the mean. In the weak coupling regime, it shows a platform, which is consistent to 
the counterpart within the Redfield scheme, independent of the coupling strength pj|. Interestingly, it shows almost 
linear behavior by enlarging the coupling strength, shown at Fig. [5fb). This demonstrates that the power noise is 
more sensitive than the energy flux in strong coupling regime as a measuring indicator. Moreover, we study the 
finite-size scaling of the optimal noise power at Fig. [5fc), which shows the linear feature. Hence, the noise signal can 
be amplified by increasing the qubits system size. Then, we analyze the relation of optimal coupling strength with 
the system size, given by Fig. [5fc). The universal scaling feature is unraveled as <a op £Oc7V _7 , 7 = 2.00T0.10. Though 
the higher orders (e.g., skewness) are not shown, the finite-size scaling effects of the optimal value and the optimal 
coupling strength are the same as those in the noise power. As a result, we conclude that both the linear scaling 
behavior of the optimal flux fluctuations and the power-law scaling behavior of the optimal coupling strength may be 
universal for the collective-qubit system in strong coupling regime. 
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FIG. 5: (Color online) (a) Steady state behaviors of noise power with different system sizes; (b) ratio of the noise power to 
photonic energy flux; (c) finite-size scaling of the optimal noise power; (d) finite-size scaling of the optimal coupling strength, 
corresponding to the optimal noise power. The parameters are given by eo = 0, lj c = 10A, Ts = 4A and Td — 2A. 


IV. CONCLUSION 

In summary, the nonequilibrium NIB A approach is applied to study the steady state quantum transport of collective- 
qubit system in strong photon-qubit coupling regime. The collective term of indirect qubit-qubit interaction, emerging 
from the cooperative contribution of many photon excitations, is essential to establish the novel steady state distribu¬ 
tion, which is apparently deviated from the counterpart within Redfield scheme. The turnover behavior of photonic 
energy flux and high order cumulants (e.g., noise power) is analyzed for finite-size qubits system, which implies the 
optimal control of the collective quantum transport. Moreover, the linear scaling feature of both optimal energy flux 
and high order cumulants is described, even in high temperature bias limit. High order cumulants show more sensitive 
signatures than the photonic energy flux by tuning coupling strength into strong regime. Interestingly, the inversed 
power-law scaling feature of the optimal coupling strength is unraveled, which is suggested to be universal for the 
nonequilibrium collective-qubit system. It is known that to detect nonmonotonic behavior of spin-boson system is still 
a challenging problem [23|, [3C|. Here, with such scaling effect, we believe it may provide an alternative road to mea¬ 
suring such nonmonotonic behavior of spin-boson system in comparatively weak coupling regime, by only increasing 
the qubits number. 
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Appendix: Quantum master equation with full counting statistics 

The general expression of second-order quantum master equation under Born-Markov approximation is given by 

(14) 

where p} is the modified reduced system density matrix, by tracing the photonic bath degrees off the total density 
matrix, and [A x , B x \ x = A X B X — B X A_ X . Starting from the reorganized system H s = eo J z — at Eq. (|2]), 
and the modified system-bath interaction V s b(x) = + e tBx J-) with the collective photonic momentum 

B x = iJ2 kv 2X * k v ( K e LUkx5v ^ D ^ 2 b\ v — iL.c.), the population dynamics under angular basis {|j, ra)} with the counting 
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field parameter is given by 

7 TDX a r oo 

= ~ { J )2 1 dT {{ C ^ eiAm - lT 9i-iP^+C(r)e- iA ^g+P^} 
+[C*(T)e- iA ™-^g+_ lP x + C*(T)e iA ^g+P*} 

~[C*(T + x) eiAm ~ lT 9m-lPm-l +C*(T + X) e ~ iArnT 9mPm+l\ 
-[C(r - x) e ~ iAm ~ 1T 9m-iPm-i + C(T - xy A ™ T g+P* +1 }} 

where the time domain correlation function is 

C(t) = = exp{— y~' ^*‘ k ’ v ) 2 [(2 n k + 1)(1 - cos w fe r) + i sinw fe r]}, 


(15) 


the angular coefficient is g+ n = j(j + l) — m(m + l) with j = N/ 2, and energy gap is A m = .E m+ i — E m = eo — (2m-bl)£ 
with eigenvalue E m = (eo m — £m 2 ). By defining the transition rates 

A f‘°° 

«m(x) = (y Am J dTC (. T ~ x)e TiAmT , 

the equation of motion for population dynamics is simplified as 

dP x 

-j^- = “fam-1 + ^m)^m + + ^m(x)^m+l> ( 16 ) 

with = ^m(x = 0). In absence of the counting field, the modified equation of motion is returned to the standard 
quantum kinetic equation for the system populations as P x=0 = P m . 



